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5 ' ABSTRACT 

q ■ The ( 2 P 3 / 2 — » 2 Pi/2) transition of singly-ionized carbon, [C n], is the primary 

coolant of diffuse interstellar gas. We describe observations of [C ill emission 
towards nine high Galactic latitude translucent molecular clouds, made with 
the long wavelength spectrometer on board the Infrared Space Observatory. To 
understand the role of dust grains in processing the interstellar radiation field 
(ISRF) and heating the gas, we compare the [C n] integrated intensity with the 
far-infrared (FIR) integrated surface brightness for the 101 sampled lines of sight. 
We find that [C n] is linearly correlated with FIR, and the average ratio is equal to 
that measured with the COBE satellite for all high-latitude Milky Way gas. There 
is a significant decrease that was not detected with COBE in [C n] emissivity at 
high values of FIR. Our sample splits naturally into two populations depending 
on the 60/im/100/iHi surface brightness ratio, or color: "warm" positions with 
60/100 > 0.16, and "cold" positions with 60/100 < 0.16. A transition from 
sources with warm to those with cold 60/100 colors coincides approximately 
with the transition from constant to decreasing [C n] emissivity. We model the 
[C n] and far-infrared emission under conditions of thermal equilibrium, using 
the simplifying assumptions that, in all regions heated by the ISRF, the most 
important source of gas heating is the photoelectric effect on grains and the most 
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important source of gas cooling is [C n] emission. The model matches the data 
well, provided the ISRF incident flux is Xo ~ 1-6 (in units of the nominal value 
near the Sun), and the photoelectric heating efficiency is e ps 4.3%. There are 
no statistically significant differences in the derived values of Xo an d e for warm 
and cold sources. The observed variations in the [C n] emissivity and the 60/100 
colors can be understood entirely in terms of the attenuation and softening of 
the ISRF by translucent clouds, not changes in dust properties. 

Subject headings: ISM: clouds — ISM: atoms — infrared: ISM: continuum — in- 
frared: ISM: lines and bands — dust, extinction 



1. Introduction 

Diffuse interstellar gas clouds are heated by the absorption of starlight and cooled by 
the emission of spectral line radiation at far-infrared (FIR) and submillimeter wavelengths. 
Energy can be transferred from the interstellar radiation field (ISRF) to the gas by the 
photoelectric effect on dust grains (Watson 1972; de Jong 1977; Draine 1978; Bakes & 
Tielens 1994), whereby far ultraviolet (FUV) (6eV < hv < 13.6 eV) photons ionize grains 
and the ejected electrons heat the gas through collisions. In diffuse H i regions of low kinetic 
temperature (T k ~ 50 — 200 K), the "cold neutral medium" (CNM), this process is thought 
to be the dominant gas heating mechanism (see Wolfire et al. 1995). 

The most abundant form of carbon gas in the CNM is the singly-ionized carbon atom, 
C ii, because the first ionization potential of carbon is 11.3 eV and the CNM is permeated 
by FUV photons. The 158/xm ( 2 P 3/2 — > 2 Pi/ 2 ) emission line transition of [C n] should be the 
primary coolant of the CNM (Wolfire et al. 1995). There are two main reasons for this. 
First, the C n ion is relatively easy to excite by collisions given average CNM temperatures 
(T ~ 80K;Kulkarni & Heiles 1987). After the hyperfine states of hydrogen, the 2 P 3/2 
fine structure state of C n is the first excited state (hu/k = 91 K) of all gas-phase CNM 
constituents (the next state is the 3 Pi level of the O i atom, which is 228 K above ground). 
Second, once excited, the energy loss rate, given by the product hv A u \ x N CII [ 2 P3/2], is much 
higher for C n than for other CNM species. Here A u \ is the Einstein spontaneous emission 



1 Based on observations with ISO, an ESA project with instruments funded by ESA Member States (espe- 
cially the PI countries: France, Germany the Netherlands and the United Kingdom) with the participation 
of ISAS and NASA. 
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coefficient for the 2 P 3/2 — > 2 Pi /2 transition and N cn [ 2 P 3 / 2 ] is the column density of C n in the 
excited state . 

Since the FUV radiation field that ionizes dust grains (and thus heats the gas) can 
also heat the grains, one expects a connection between the [C n] 158/xm line emission from 
photoelectron-heated gas and the far-infrared (FIR) emission from FUV-heated dust grains. 
Large-scale surveys of the Galaxy indeed reveal a link between the C n and FIR emission 
properties of the CNM. An unbiased survey of spectral line and continuum emission in the 
Milky Way was conducted with two instruments on board the Cosmic Background Explorer 
(COBE) satellite: the Far Infrared Absolute Spectrometer (FIRAS) and the Diffuse Infrared 
Background Experiment (DIRBE). The [C n] emission was observed with FIRAS to be closely 
correlated with H i emission at high Galactic latitudes, with an average C n cooling rate of 
(2.65 ± 0.15) x KT^ergs- 1 (H atom)" 1 (Bennett et al. 1994). Both DIRBE and FIRAS 
maps of the dust continuum show that the dust opacity and FIR surface brightness are also 
correlated with H i column density at high latitude (Boulanger et al. 1996). 

So on the large scale (FIRAS has an angular resolution of 7°, DIRBE of 40'), CNM 
H i gas has constant FIR and [C n] emissivity. It is therefore likely that the processing 
of FUV radiation by dust grains (1) to produce FIR dust emission and (2) to heat the 
gas — resulting in [C n] emission — occurs in similar ways throughout the cold interstellar 
medium. It is difficult to test this idea directly using the COBE data, mainly because of 
inadequate resolution. The emission from gas with vastly different physical conditions, for 
example molecular and atomic gas, might be combined together in a single FIRAS or DIRBE 
measurement, making it impossible to understand the physical mechanisms responsible for 
the observed emission. 

High resolution observations of the [C n] and FIR emission from translucent molecular 
clouds can be used to test the limits of models for gas heating by the ISRF. Molecular gas 
can only exist in shielded locales where the FUV portion of the ISRF has been attenuated, 
and thus where the photoelectric heating is much weaker than in atomic regions. In these 
regions, most of the carbon has combined into CO, and most of the [C n] emission has 
"turned off." Although FUV radiation is weak in molecular zones, optical photons are 
still present to heat dust grains and produce FIR emission. Indeed, many IRAS cirrus 
clouds (Low et al. 1984) first detected in FIR emission at 60 and 100/im are associated 
with intermediate-extinction molecular clouds, or translucent clouds (van Dishoeck, & Black 
1988). Translucent molecular clouds are most easily observed when nearby and at high 
Galactic latitudes {b > 15°) (Magnani, Blitz, & Mundy 1985). Currently over 100 nearby 
((d) ps 105 pc) translucent high-latitude clouds (HLCs), have been cataloged using their CO 
emission (see Magnani, Hartmann, & Speck 1996). In a study of 75 HLCs, it was shown 
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that all of the clouds are associated with H i, and that the HLCs probably condensed out of 
diffuse CNM gas (Gir, Blitz, & Magnani 1994). 

We examine here the [C n] and FIR emission towards a sample of translucent HLCs 
using the ISO satellite. Since the mean angular size of HLCs in CO emission maps is ~ 1° 
(Magnani, Hartmann, & Speck 1996), the clouds are easily resolved by the ~ 71" beam of 
ISO. This allows us to choose between lines of sight towards individual clouds with varying 
amounts of atomic and molecular gas, which was not possible with COBE. As described 
above, the lines of sight dominated by molecular gas include regions from which FUV photons 
are excluded. 

In this paper we summarize the ISO observations of [C n] emission towards HLCs, all 
of which are available from the ISO archive (§2). We discover that HLCs fall into two 
categories based on their 60/xm/100/im color, "cold" sources and "warm" sources. We find 
that the line-of-sight [C n] emissivity is consistent with the average high-latitude COBE 
cloud emission for both cold and warm positions. We surmise that the chief distinction 
between the two 60/100 color regimes is the column density (§3). Assuming that the local 
heating and cooling depend only on the fixed properties of dust grains and the attenuated 
ISRF intensity, we develop a model for the [C n] and FIR emission of translucent clouds. 
The model enables us to estimate the average line of sight photoelectric heating efficiency, 
e, as well as the intensity of the interstellar radiation field incident on HLCs, xo (§4). We 
discuss the significance of our results, and point out some issues unresolved by this study 
(§5). Finally, we summarize our conclusions (§6). 



2. Observations 

2.1. ISO Observations of Translucent High-Latitude Clouds 

The Infrared Space Observatory (ISO) (Kessler et al. 1996), a satellite observatory 
built by the European Space Agency (ESA), was launched in 1995 November and operated 
until its cryogenic reserves were depleted in 1998 April. A total of seven peer-reviewed 
observing programs dedicated to studying [C n] emission from nine translucent high-latitude 
molecular clouds were carried out with the long wavelength spectrometer (LWS) (Clegg et al. 
1996) on board ISO. Some of these [C n] measurements were the results of proposals by us 
(observer:proposal = TBANIA:TMB_1, TBANIA : TMB_2, TBANIA:SH_HLC, and WREACH : PHASE) , 
and the rest were extracted from the ISO public data archive. 2 The TBANIA data were 



2 See the web site http://www.iso.vilspa.esa.es/ for details. 
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originally published in Ingalls (1999). The data for twenty of the positions observed towards 
cloud MBM-12 (L1457) were reported previously by Timmermann, Koster, & Stutzki (1998). 

Table 1 lists the nine HLCs observed in [C n] with ISO, together with their positions 
in Galactic coordinates (£,b) (Columns 2 and 3). Each entry in the Table represents a 
separate source "observation," i.e., a self-contained set of measurements of the target. A 
given observation may consist of a single point, or it may consist of multiple pointings of 
the telescope as part of a raster map. For clouds with \b\ > 25° (all clouds in our sample 
but G300. 1-16.6) we use the cloud names given in the Magnani, Hartmann, & Speck (1996) 
catalog, and list commonly used alternate names in parenthesis (column 1 of Table 1). The 
field of view (FOV) of each observation (column 4) is defined by the full width at half 
maximum (FWHM) beam size of the ISO LWS instrument as well as the total extent of 
the raster grid employed. We indicate the position angle (PA) of raster grids, measured 
counterclockwise from celestial North, in column 5. 

Emission spectra of [C n] were measured with the LWS spectrometer in either LWS01 
or LWS02 scanned grating modes. Rapid scanning was used, and each spectral resolution 
element (of width ~ 0.6yum) was sampled many times. This ensured that cosmic ray hits 
on the detector were efficiently diagnosed and removed. We used detector LW4, centered 
on 161 /im, to measure the 157.7/zm ( 2 P 3 / 2 — ► 2 Pi/2) transition of [C n]. This detector is 
estimated to have a nearly circular beam, with FWHM size ~ 71" (Swinyard et al. 1998), 
20% smaller than predicted by early optical models. The truncation of the LWS beam is 
attributed to multiple clipping effects and stray light propagation within the optical system 
(Caldwell et al. 1998). Stray light is also apparently the cause of interference fringes, which 
manifest as a sinusoidal modulation with period ~ 10/im at A = 158/xm in the LWS spectra 
of extended sources (Trams et al. 1996). We attempted to minimize the effect of the fringes 
on our [C n] integrated intensities by subtracting either a linear or quadratic baseline from 
the spectral data measured with the LW4 detector. 

We derived C n integrated intensities using the Automatic Analysis (AA) data products 
from the ISO archive. All calibration steps were performed at ESA in batch mode, i.e., 
without human input. In order to check the quality of the AA procedure, we calibrated a 
subset of our spectra manually using the LWS Interactive Analysis software and compared 
the results with the the corresponding AA results. The integrated [C n] intensities produced 
using the two methods were indistinguishable to within ±10%. As the absolute flux calibra- 
tion of the LWS has been estimated to be no better than ±20% (Burgdorf et al. 1998), we 
judged that calibration beyond the AA stage was not necessary. 
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We reduced the calibrated scan data using the ISO Spectral Analysis Package (ISAP). 3 
First, we averaged all scans and scan directions for detector LW4 in bins of width 0.2/zm using 
a 3o" median clipping algorithm. Then we subtracted linear or quadratic baselines from the 
continuum, and made Gaussian fits to each spectrum. Since the Doppler equivalent velocity 
resolution of the spectra is ~ 1000 km s -1 , which is much larger than the typical H i linewidth 
of HLCs (6-26 kms" 1 ), we assumed that all measurements were unresolved spectrally. This 
allowed us to fix the FWHM linewidth in the Gaussian fits to the known LWS resolution of 
0.6/iin. We used the Gaussian fits to compute the [C n] integrated flux, which we divided 
by the LWS beam area, O = 9.3 x 10~ 8 sr, to obtain the [C n] integrated intensity, I c „, in 
units of erg cm -2 s" 1 sr" 1 . We derived error bars in I CII by adding a 20% absolute calibration 
uncertainty (Burgdorf et al. 1998) in quadrature with the formal Gaussian fit errors in the 
[C n] line integral, obtained using the ISAP software. 



2.2. IRAS 60 and lOO^m Measurements 

We obtained IRAS Sky Survey Atlas (ISSA) images of each cloud and extracted the 
60/iin and 100/xm surface brightness, I 60 and I 100 respectively, for the positions observed in 
[C n]. The angular resolution of the ISSA is only 4', whereas the ISO LWS has a FWHM 
beamsize of 71" at 157.7/im, so we interpolated the ISSA maps onto a 71" grid before making 
I 60 and I 100 measurements. 

We have tested the validity of interpolating the 4' IRAS images to compare with ISO 
data. We obtained HIRES-processed IRAS images of the clouds in our sample (Levine et 
al. 1993). The HIRES maps have approximately 60" and 100" resolution at 60 and lOO^um, 
respectively, which is comparable to the ISO LWS beamsize. A linear fit to the HIRES data 
as a function of interpolated ISSA surface brightness shows a strong correlation with slope=l 
for pixels with ISSA values brighter than about half of the maximum value. Below the half- 
maximum threshold, HIRES images have an average surface brightness of zero, probably 
the result of background subtraction in the HIRES process. We conclude that on average 
the interpolated ISSA images give a good approximation to the actual surface brightness 
at higher resolution, with the advantage of maintaining the low surface brightness emission 
that does not appear in HIRES images. 

The ISSA plates are optimized for relative, not absolute, photometry (Wheelock et al. 
1994). To put images of different clouds on the same relative scale, we devised a zeropoint 



3 The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument 
Teams and Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, RAL and SRON. 
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adjustment method that took advantage of the constant 60/im/lOO/iin scaling properties 
of high-latitude clouds. Pixel-by-pixel comparisons of the 60/mi and 100/rni images show 
a strong correlation (the correlation coefficient is typically greater than 0.9). Linear fits 
to the I 100 versus I 60 data for each of our nine clouds yield the average relationship I 100 = 
(5.74 ± 0.15)I 60 - (1.14 ±0.13). The relative variation in the derived slopes is 0.15/5.74~ 3%, 
showing that the 60/100 ratio is nearly constant for all cloud images. On the other hand, 
the relative variation in the I 100 -intercepts exceeds 11%, indicating significant fluctuations in 
the 60/im and 100/im surface brightness zeropoint from image to image. We adjusted the 
zeropoint of each image so that I 1( ,„=0 when I 60 =0 by subtracting the derived I 100 -intercept 
values from each 100/im cloud image. 

We converted the IRAS relative scale to the COBE DIRBE absolute surface brightness 
scale by multiplying the 60/im data by 0.87 and the 100/im data by 0.72 (DIRBE 1998). 
We derived error bars in I 100 by combining in quadrature the 10% accuracy in point source 
flux recovery (IRAS 1988), with the 7% uncertainty in the DIRBE to IRAS conversion at 
100/im. We derived error bars for I 60 measurements by combining in quadrature the 10% 
accuracy in point source flux recovery with the 2.5% uncertainty in the DIRBE to IRAS 
conversion at 60/im. 

3. Observational Results 

A total of 109 positions in HLCs were observed using the ISO LWS, and [C n] emission 
was detected towards 101 of them. We list the [C n] detection rate for each observation in 
column 6 of Table 1. 



3.1. The HLCs as Representative [C n] Clouds 

Ingalls et al. (2000) showed that HLCs are the average Galactic high-latitude sources 
of CO and far-infrared emission. We demonstrate here that they are also representative 
sources of [C n] radiation. We mentioned in §1 that CNM H i emission is correlated with 
both far-infrared and [C n] emission. Now we compare directly the FIR and [C n] properties 
of high-latitude material using the all sky COBE surveys. The FIRAS Line Emission Maps 
give [C n] 157.7/im integrated intensity at 7° spatial resolution. The DIRBE Zodi-Subtracted 
Mission Average (ZSMA) maps give dust continuum surface brightness in far-infrared bands 
at 40' resolution (DIRBE 1998). We use a function introduced by Helou, Soifer, & Rowan- 
Robinson (1985) (see also Helou et al. 1988) to define the far-infrared surface brightness 
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in terms of the 60 and 100/zm surface brightness: 



FIR = U 6 A„ + 2.58U . 



10~ 6 erg cm -2 s" 1 sr 1 \ MJy sr 

This representation of far-infrared emission is estimated to be accurate to within 1% for 
modified blackbody spectra from dust clouds with temperatures between 20 and 80 K and 
emissivity that varies as z/ " 2 . 

We created DIRBE maps of FIR using the 60 and 100/im ZSMA data, resampled to 
the 7° FIRAS grid, and compared with the [C n] map (see Figure 1). We obtain a highly 
significant correlation (there is less than 0.05% chance that the data are uncorrelated) . A 
linear fit weighted by the error bars for all data with |6| > 5° gives 

I CII = (2.54 ±0.03) x 10" 2 FLR- (1.1375 ±0.0004) x 10" 6 erg cm" 2 s" 1 sr" 1 , (2) 

where I CII and FIR are in the same units. A grey wedge on Figure 1 shows the 3a 
range of expected FIR and I CII values. Since the fit was weighted by the COBE error 
bars, this wedge is much narrower than the actual scatter in data points (the rms scat- 
ter in I CII about the fit is 5.29 x 10~ 6 ergcm~ 2 s" 1 sr" 1 ). The fit gave an I CII -intercept of 
— 1.1375 x 10" 6 erg cm" 2 s" 1 sr -1 , because many of the FIRAS [C n] values for low FIR are 
negative. For what follows, we assume that in the absence of FIR emission there is no [C n] 
emission, and ignore the intercept. 

We compare the [C u]-FIR relationship for "COBE clouds" with our HLC data. The 
C ii integrated intensities of the 101 HLC positions detected with ISO are plotted as a 
function of the FIR surface brightness in Figure 2. Linear regression to the data, weighted 
by the errors, gives 

Ion = (2.5 ± 0.9) x 10~ 2 FIR. (3) 

In the fit we have fixed the I CII -intercept at zero to reflect our assumption that I CII = 
when FIR = 0. We plot the fit as a solid line in Figure 2. A gray wedge labeled "COBE" 
represents the ±3<r range of expected (FIR, I CII ) values based on the COBE data (Equation 
2). The HLC and COBE slopes are identical. Our analysis is consistent with the HLCs 
being part of the same population of sources as the average CNM source observed with the 
low-resolution COBE instruments. 

There is evidence that at the ISO resolution a straight line is not the best representation 
of the l CII -FIR data. Breaking the data into FIR subregions and performing separate linear 
fits yields a systematic decrease in slope as FIR increases. Figure 3 shows an example of 
this for bins of width 10 -4 erg cm -2 s _1 sr" 1 along the FIR-axis. The local slope is close 
to the overall slope for FIR < 2 x 10~ 4 erg cm -2 s" 1 sr" 1 , but is virtually zero for FIR > 
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2 x 10~ 4 ergcm^ 2 s -1 sr -1 . We also plot in Figure 3 the weighted mean values of I CII and 
FIR in the bins. A significant departure from the COBE and ISO predictions is seen in 
the bin with the largest FIR integrated surface brightness, where I CII is 25% lower than 
expected. Thus, the [C u] emissivity of HLC gas decreases as FIR increases. In §4 we 
explain this using a model of the radiative heating and cooling of translucent clouds. 

3.2. The 60/100 color 

Ingalls et al. (2000) showed that the eight HLCs they studied have typical 60/100 colors 
(I 60 /I 100 ) when compared to COBE clouds. The 60/100 colors for our nine clouds are likewise 
typical. We plot I 100 as a function of I 60 for our sample in Figure 4. Superimposed on the 
plot is a dashed line representing the average behavior for Milky Way gas correlated with H i 
emission, I 60 = 0.16 1 100 (Dwek et al. 1997). Most positions fall slightly below this line. These 
sources probably contain mostly atomic gas, since their FIR surface brightnesses are low 
and their I 60 /Ii o ratios are slightly "warmer" than that of average H i. For I 100 > 4MJysr _1 
there is a split into two populations: the "cold" sources with I 60 < 0.16 1 100 ; and the "warm" 
sources with I 60 > 0.16I 1( ,„. 

Warm positions generally produce less emission than cold positions — over 90% of the 
warm positions, but less than 30% of the cold positions, have FIR < 2xl0~ 4 ergcm~ 2 s" 1 sr _1 . 
Their [C n] and FIR properties are not easily distinguishable otherwise. Weighted linear fits 
to the [C n] vs. FIR data, holding the intercept fixed at zero, give slopes of (2.6± 1.1) x 10~ 2 
for warm sources and (2.2 ± 1.7) x 10~ 2 for cold sources. These slopes are equivalent within 
the errors to each other, and to the slope determined for the sample as a whole. 

The most striking difference between cold and warm positions seems to be their column 
density. Figure 5 shows this for high-latitude cloud MBM-12 (L1457). Here we compare 
the 60 and 100/xm images to a visual extinction map of the cloud. The extinction, A v , 
was derived using an adaptive grid star count algorithm (Cambresy et al. 1997) from H 
and K images in the 2MASS second incremental data release. 4 After resampling the 60/im, 
100/im, and A v images onto the same grid, we compared the 60/100 color with A v for each 
pixel. Figure 5 shows I 60 /Iioo averaged in bins of width AA V = 0.71 mag. The 60/100 color 
decreases monotonically with increasing extinction, implying that warm dust emission is a 



4 This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project 
of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute 
of Technology funded by the National Aeronautics and Space Administration and the National Science 
Foundation. 
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characteristic of low extinction positions, while cold dust emission is a characteristic of high 
extinction positions. The effect is not slight. A small adjustment in 60/100, e.g., from 0.184 
to 0.160, represents a large change in A v , from approximately 1 to 3. 

In what follows, we develop a model for the radiative transfer of the ISRF in translucent 
clouds. This allows us to predict the emission from dust grains, as well as the photoelectric 
heating and [C n] cooling of the gas. We show that the cold/warm dichotomy in translucent 
clouds is probably caused by the attenuation of the ISRF, not changes in dust properties. 

4. Modeling the C n and FIR emission of HLCs 

The [C n] emissivity towards translucent HLCs decreases as the FIR surface brightness 
increases. The cold and warm populations of HLC positions, identified via their 60/100 col- 
ors, have different column densities, and hence different amounts of atomic versus molecular 
gas, but similar [C n]-FIR emission properties. To explain these empirical facts we offer a 
model of the heating and cooling of gas and dust in the interstellar medium. The model 
describes regions heated primarily by the interstellar radiation field, i.e., translucent regions. 
We simplify the calculations by treating only the single most important sources of heating 
and cooling of translucent gas: the photoelectric effect for gas heating, and [C n] emission for 
gas cooling (Wolfire et al. 1995). We predict the FIR and [C n] emission under conditions 
of thermal equilibrium. For the dust grains, this requires that the power emitted by each 
grain equals the power absorbed from the attenuated interstellar radiation field. For the gas 
this requires that the local [C n] cooling rate equals the photoelectric heating rate. 

By assuming that [C n] emission is tied to photoelectric heating, we bypass chemistry. 
We do not attempt to determine the abundance of C n gas in the [C ii]-emitting regions, 
e.g., by balancing the production and destruction of C n using a chemical model. Instead 
we determine the photoelectric heating, which depends on the dust properties, and equate it 
with the [C n] cooling modified by an efficiency parameter (e.g., see Bakes & Tielens 1994). 
This approximation assumes that (1) most [C n] cooling takes place in FUV-heated areas 
and (2) most FUV heating takes place in areas cooled by [C n] emission. 

If C ii were abundant in regions not heated by the photoelectric effect, then much of the 
[C n] emission might be the result of another heating mechanism, such as cosmic rays, and 
we would underestimate I CII . This is unlikely. If the radiation field is too weak to ionize dust 
grains (average work function: 5.5 eV), it is also too weak to ionize C i (ionization potential: 
11.3 eV) and form C n. In other words, where there is [C n/ cooling, it is the result of FUV 
heating. 
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The second tier of our approximation maintains that most photoelectric heating occurs 
in regions cooled by [C n] emission. If this did not hold, then the emission from another 
species like C i or CO might provide most of the cooling in some FUV-heated regions and our 
I CII predictions would be overestimates. It is certainly true that C i and CO are important 
coolants of translucent clouds, but in terms of total flux the [C i] and CO lines are much 
weaker than the [C n] line. The brightest known HLC source of [C i], MBM-12, has a 
maximum observed [C i] ( 3 Pj — > 3 P ) integrated intensity of about 3.6 x 1CT 7 erg cm~ 2 s _1 sr _1 
(Ingalls, Bania, & Jackson 1994), a factor of 10 weaker than the [C n] emission towards most 
of our HLCs. According to the photodissociation region models of Kaufman et al. (1999), 
the other [C i] fine structure transition, ( 3 P 2 ^ 3 -Pi), should have at most twice the power 
of the ( 3 P! — > 3 P ) transition under translucent conditions (FUV flux equal to the local 
interstellar value and gas volume density ~ 10 2 ~ 5 cm~ 3 ). Even if all of the C i atoms that 
emit in these transitions were excited by FUV-ejected electrons, the photoelectric heating 
in the [C i]-emitting regions of HLCs is still at most only 30% of that in the [C ii]-emitting 
areas. Furthermore, CO cooling is almost an order of magnitude weaker than [C i] cooling 
for translucent clouds (eg., see Ingalls et al. 2000). Thus, where there is significant FUV 
heating, most of the cooling is by [C n] radiation. 

4.1. The penetration of the interstellar radiation field into translucent clouds 

Here we calculate the intensity spectrum of the interstellar radiation field inside a model 
cloud, as it is attenuated and processed by dust grains. We adopt a plane parallel geometry, 
where changes in the mean radiative properties of the gas and dust are dependent only on 
the optical depth, r A , measured from the surface of the cloud to the center. The cloud size 
is defined by the central optical depth, r c> \. 

The interstellar radiation field, Ji SR f,a, is isotropic and incident on both sides of our 
plane parallel clouds. We consider only the stellar component of the ISRF longward of 
0.0912/im (the Lyman continuum), and the cosmic background radiation field. We use the 
ISRF approximation of Mezger, Mathis, & Panagia (1982) and Mathis, Mezger, & Panagia 
(1983), which is the sum of a stellar UV component and three diluted blackbodies. We add 
to this the cosmic background radiation field, _Ba(2.7K), where B\(T) is the Planck function. 
We calculate the integrated intensity of the interstellar radiation field to be 



*/0.0912/zm 

This is 27% higher than the value quoted by Mathis, Mezger, & Panagia (1983), partly 
because we consider a larger range of wavelength than they did. In constructing an ISRF, 




(4) 
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we do not include clouds themselves as infrared sources. Infrared photons have little affect on 
the energetics of translucent clouds because the clouds are transparent to them. Furthermore, 
the interstellar infrared spectrum results from the processing of stellar radiation by clouds, 
and does not represent a new source of energy. A goal of this paper is to calculate the 
infrared spectrum of translucent clouds, to compare with observations. 

We use the technique of Flannery, Roberge, & Rybicki (1980) to model the penetration 
of the ISRF into interstellar clouds. They expressed analytically the solution of the radiative 
transfer equation in the presence of scattering using a spherical harmonics method. The 
mean intensity of the radiation field at depth r A , for a cloud of size t CjA is: 

M 



m=l 



x {exp[-/c miA t x ] + exp[-/c mjA (2r CjA - r A )]} . (5) 

The units of J A are ergs -1 cm -2 sr _1 /jm" 1 . The coefficients, A m> x{xo), are determined using 
the boundary condition at r A = 0, i.e., J A (0, r CjA , xo) = Xo <Asrf,a, whence: 

[1 + exp(-2fc m , A T C! x)\ M 

The multiplicative factor, x 0! allows for variations in the nominal radiation field (Mezger, 
Mathis, & Panagia 1982; Mathis, Mezger, & Panagia 1983). The scattering modifiers, 
{k m ,x}, are the reciprocal eigenvalues of an M x M tridiagonal symmetric matrix, X, given 
by Equation (A3) of Flannery, Roberge, & Rybicki (1980). If we use the scattering phase 
function introduced by Henyey & Greenstein (1941), then the matrix elements are functions 
only of the dust grain albedo, (u)\), and the mean cosine of the scattering angle, (g\). The 
() notation we use denotes an average of dust optical properties over the grain population 
(see Appendix). The Ith off-diagonal element of the matrix, X/^+i, equals // \fhi~Jii, where 

h l = (2l + l)(l-(u x )(g x ) 1 ). (7) 

The matrix is symmetric, so A; +lj ; = A^ +1 . In practice we truncate the series in Equation 
5 to M = 20 terms, since for the wavelengths of interest foi is less than 1% of the value of 
k 2 o- Therefore, X is a 20 x 20 matrix. 

Equation 5 takes as inputs the optical depth, r A , and the total size of the cloud, given 
by t c A , but we can switch our depth scale to that of the visual extinction, A v = 1.086r v . The 
visual extinction is a convenient depth scale because it can be related linearly to the column 
density of hydrogen nuclei, N H (eg., Bohlin, Savage, & Drake 1978). It can be defined in 
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terms of the wavelength-dependent optical depth using the extinction curve (see Equation 
A2): 

A v = 1.0867-a (J^J . (8) 

The central visual extinction, A V)C , can be defined similarly in terms of t c< \. 

We have used Equation 5 to calculate the mean intensity spectrum inside interstellar 
clouds (Figure 6). In what follows, we use the Equation 5 attenuated radiation field to 
compute the heating function for dust grains and the photoelectric heating function for the 
gas. Our assumption of thermal equilibrium then allows us to estimate the surface brightness 
of 100/im and 60/im radiation from the dust and the [C n] emission line cooling of the gas. 



4.2. Dust and gas heating 

4-2.1. Thermal equilibrium of dust grains 

The temperature of a spherical grain of radius a and material % (see Appendix), measured 
at depth A v in a cloud of central thickness A v c , is T g (a, i, A v , A v c , xo)- It can be determined 
by equating the power absorbed by the grain to the power emitted. This requires: 

poo poo 

4ttV / J x Q abs , x (a,i)d\ = 4ir 2 a 2 B x (T g ) Q abs , x (a,i) dX, (9) 
Jo Jo 

where B X {T) is the Planck function and Q a bs,\(a,i) is the absorption efficiency (see Ap- 
pendix). The assumption here is that dust grains are in thermal equilibrium. This holds 
for grains with sizes a > 250A, but breaks down for smaller particles, which are stochasti- 
cally heated (See Draine & Li 2001, and references therein). We show in §5.3.1 that this 
equilibrium approach does not affect significantly our conclusions. 

We have solved Equation 9 numerically for the function T g (a, i, A v , A v c , xo)- As an ex- 
ample, we plot in Figure 7 the equilibrium temperature of a graphite (i = 1) grain of size 
a = 0.1/iin as a function of A v , for a cloud with A ViC = 1.0 mag, immersed in a xo = 1 inter- 
stellar radiation field. At the surface of the cloud, a 0.1/zm graphite grain has an equilibrium 
temperature of T g m 21 K. This is comparable to the value independently computed by Li & 
Draine (2001) for "diffuse clouds," i.e., clouds without ISRF attenuation. For all grain sizes 
> 0.01/iin and for both graphite and silicate particles, our values of T g at the surfaces of 
clouds match those of Li & Draine . Our model, which computes the grain temperatures as 
a function of cloud size, depth, and xo-, can be used to extend the Li & Draine model to the 
translucent regime, where radiative transfer is important. Figure 6 shows how the spectral 
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shape of the radiation field is already significantly different at the center of an A VjC = f cloud 
from that at the surface. Obviously for translucent clouds, a depth-dependent approach to 
computing grain temperatures is necessary. 

4-2.2. Photoelectric heating of gas 

The most important way to transfer energy from the ISRF to the gas is via the photoelec- 
tric effect on dust grains. Photoelectric heating of the gas can occur when photons with ener- 
gies larger than the grain work function expel electrons. A typical neutral grain work function 
is ~ 5.5 eV, so it is mainly FUV photons (6eV < hv < f3.6eV ; 0.2066//m > A > 0.0912/xm) 
that contribute to the photoelectric effect. The photoelectric heating mechanism has been 
explored in detail by Bakes & Tielens (1994). They found that a fraction, e, of the total 
FUV energy absorbed by all grains is available as photoelectron kinetic energy to heat the 
gas. The value of e, the "photoelectric heating efficiency," is expected to be ~ 3% for neutral 
grains (Bakes & Tielens 1994). This can be explained using a simple heuristic argument. 
An average FUV photon has energy of 8 eV, which causes an electron with 2.5 eV of kinetic 
energy to be released from a neutral grain with work function 5.5 eV. In other words, 30% of 
the incoming photon's energy is liberated. The typical photoelectric yield is 10%, i.e., 10% 
of all photons absorbed by grains actually ionize them. Therefore, the fraction of photon 
energy absorbed that gets injected into the gas as kinetic energy is e ~ 0.1 x 0.3 = 0.03. 

The neutral grain photoelectric efficiency is an upper limit to the actual value of e. 
For ionized grains, the work functions are higher than 5.5 eV, and e is less than the neutral 
grain limit. In the cold neutral medium, less than 50% of grains with a < 15A are expected 
to be ionized (Li & Draine 2001; Weingartner & Draine 2001), implying relatively high 
photoelectric efficiency in translucent clouds. 

We used our cloud model to estimate the gas heating in translucent clouds. To compute 
the total amount of FUV energy absorbed by grains we first estimated the mean absorption 
cross section per H atom as a function of wavelength, (cr\): 

( a *) = y2 ° 7ra 2 C i g a6s , A (a,i)a- 3 - 5 da, (10) 
i=i j3 - 2A 

where {Cj} and a max are defined in the Appendix. Bakes & Tielens (1994) estimated that 
approximately one half of the heating in interstellar clouds is caused by tiny carbon grains 
and Polycyclic Aromatic Hydrocarbons (PAHs) with sizes less than 15A, and the other half 
comes from larger grains with 15A < a < 0.01/xm. The size distribution in the Mathis, 
Rumpl, & Nordsieck (1977) model (Equation Al) has a minimum cutoff, a min = 0.005/im, 
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so we had to extend the grain parameter calculations of Draine & Lee (1984) to include 
smaller graphite grains. Following Bakes & Tielens (1994), we assumed that for grains with 
o < A (the Rayleigh limit) the absorption and scattering efficiencies for graphite grains, 
Qabs,x{ a , 1) an d Qsca,\( a , l)j respectively, scale with grain volume. We also assumed that 
the mean scattering cosine for graphite, g\(a, 1), follows the Mie expansion given in Laor & 
Draine (1993, their Equation 8). 

The local gas heating function, T, is derived by integrating the mean intensity of the local 
radiation field, weighted by the absorption cross section per H atom and the photoelectric 
efficiency, over all FUV wavelengths and solid angles: 

r = 47re / Jx(ax)d\. (11) 

^0.0912/im 

The units of T are erg s -1 (H atom)" 1 . In Figure 7 we show the behavior of T as a function 
of A v for a cloud with A v c = 1.0 mag, e = 0.043, and xo — 1- Near the surface of the cloud 
(A v ~ 0.3 mag), the heating from our model equals the average cooling observed with COBE 
for high-latitude H i gas, 2.65 x 10~ 26 ergs -1 (H atom) -1 (Bennett et al. 1994). 



4.3. Dust and gas cooling 

The grain temperature and the gas heating function allow us to derive the FIR intensity 
and the [C n] intensity, respectively, towards translucent clouds to compare with observations 
of these important cooling sources. In the context of our proposed model these are the sole 
mechanisms of cooling. 



4.3.1. FIR Emission 

Integrating over the grain distribution, the local cooling function of dust by thermal 
emission is: 

A djA = J2 B AT 9 ) ™ 2 Qa bs ,x(a, i) a" 3 - 5 da, (12) 

i = \ >' a min 

where 

Bu{T) _ 3.97 x 10- 11 

MJysr- 1 exp(1.44 x 10~ 4 /iin • K/AT) — 1 1 ' 

is the Planck function. Since B U [T] is in MJy sr -1 , the units of A d are MJy sr" 1 - cm 2 (H atom) -1 . 
This is equivalent to MJysr -1 per H column density, N H . The emergent thermal dust emis- 
sion spectrum can be derived by converting to MJy sr -1 per A v using the standard relation, 
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N H (cm 2 ) = 1.87 x 10 21 A V (mag), integrating A d from the surface of the cloud to the center, 
and multiplying by 2: 

J A (A V , C , xo) = 3.74 x 10 21 / A d , A (A v , A V)C , X o) dA v . (14) 

Jo 

The factor 2 is necessary because the integral is from the edge to the center of a symmetrical 
one-dimensional cloud. 

4.3.2. FIR Spectra 

Model FIR spectra, I\(A VjC , xo), are plotted in Figure 8 for xo — 1- To check our model 
for internal consistency, we realize that the integrated FIR emission should be less than or 
equal to the incident ISRF integrated intensity, 2.76 x 10~ 2 erg cm -2 s^ 1 (Equation 4). For 
A v>c = (0.5, 1.0,2.0,5.0), the integrated FIR flux from the clouds is (1.00, 1.53,2.03,2.64) x 
10 _2 ergcm _2 s _1 . Thus, only for extremely opaque clouds (2A VjC ~ 10.0 mag) is most of the 
input radiation recovered as FIR emission, and the input flux is never exceeded. 

We can use Figure 8 to predict that clouds with large A v c appear "colder" than clouds 
with small A v c , provided xo is held constant. The long wavelength spectrum mimics the 
emission by dust at a single temperature. Thick gray curves in Figure 8 show modified 
blackbody spectra for clouds with the same sizes (values of A VjC ) as the model clouds, but 
with different "dust temperatures," T d =(19,18,16.9,15.1) K, for A ViC = (0.5,1.0,2.0,5.0). 
Here the blackbody function at Td was modified by the dust emissivity, assumed to be 
proportional to A~ 2 . The modified blackbody curves reproduce the model spectra quite well 
for wavelengths longward of the peak. 

We find that a thicker cloud produces more emission than a thinner cloud, but its spec- 
trum peaks at a longer wavelength (vertical arrows), yielding a lower "dust temperature," as 
shown by the gray curves. The same behavior was noticed for a large sample of high-latitude 
clouds by Reach, Wall, & Odegard (1998), who found that the average color temperature 
for dust mixed with HLC atomic gas is (17.0 ± 0.3) K, whereas for dust mixed with HLC 
molecular gas the average color temperature is (14.7 ± 0.2) K. This happens because, even 
though opaque clouds absorb more of the ISRF than transparent clouds, most of the dust 
in their interiors is exposed to an attenuated version of the ISRF that is deficient in high 
energy photons. 

Our model can be used to estimate the 60 and 100/xm surface brightness, I 60 and I 100 , 
towards translucent clouds, as well as the IRAS^derived FIR integrated surface brightness 
(Equation 1). Figure 9 plots I 60 , I 100 , and FIR versus A V)C for xo — 1- Note that even 
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though FIR increases with increasing A VjC , the rise is not linear, i.e., the FIR emissivity 
(F/i?/N H ) of translucent clouds decreases with increasing N H , as is observed (see Reach, 
Wall, & Odegard 1998). 

Another way to demonstrate that thicker clouds have "cooler" dust emission is to com- 
pare the surface brightness ratio I 60 /Iioo with visual extinction, as we have done empirically 
for MBM-12. Gray curves on Figure 5 show that, in our model, the 60/100 color decreases 
with A v (= 2A V)C ) and increases with xo- The model agrees remarkably well with the MBM- 
12 data for xo — 1-6 ± 0.2. Such behavior was already noticed and discussed qualitatively 
for translucent clouds in the Lynds 134 cloud complex (which contains two objects in our 
sample, MBM-37 and MBM-33) (Laureijs, Clark, & Prusti 1991; Laureijs et al. 1995). 
Cuts across the densest regions of L134 (Laureijs et al. 1995) show that the 60/100 color 
decreases as B-band extinction, A B , increases. We have compared the Laureijs et al. L134 
data to our model (using 1AA B ~ A v ), and we find a good match for xo — 1-4 — 1.6 and 
A v = — 5. There is a sudden drop in 60/100 for extremely opaque (A v > 5) lines of sight 
towards the L134 complex (Laureijs, Clark, & Prusti 1991; Laureijs et al. 1995), which the 
model does not predict. The authors attribute this to a change in dust properties (see §5.3). 



We estimate the [C n] integrated intensity towards translucent clouds. Recall our key 
assumptions that the gas cools only through emission in the ( 2 P 3 / 2 — > 2 Pi/2) transition of 
C ii, and that the gas is in thermal equilibrium. We set the cooling, A (i.e., the [C n] emission 
per N H ), equal to the gas heating, T, given by Equation 11. The integrated intensity of the 
[C n] emission is: 



We plot in Figure 9 the [C n] integrated intensity, I en, as a function of A v c for clouds 
with xo — 1 and e = 0.043. The two upper curves highlight the differences between gas and 
dust heating by interstellar radiation. The qualitative behavior of both curves is the same. 
That is, both FIR and I CII rise almost linearly, implying linear increases in column-averaged 
dust and gas heating with extinction; then both curves show a decrease in slope, implying 
that clouds with central extinction larger than a transition size produce much less emission 
at depths greater than this size. The effect of reaching the transition sizes is different for the 
two tracers. The I CII emission becomes flat at A v c ~ 0.5 mag, whereas FIR only shows a 
decrease in slope, without flattening even for A v c = 2.0 mag. In the interior of a translucent 



4-3.3. [C n] emission 




(15) 
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cloud the mean integrated intensity of FUV radiation is much smaller than that of optical 
radiation (refer to Figure 6). Dust grains are heated by ISRF photons of all wavelengths, 
but photoelectrons, which lead to [C n] emission, are ejected only by FUV photons. Thus, 
one expects the [C n] emission to extinguish before the FIR emission does. 

4-3.4- Linear and Flat Behavior 

The model for FIR and [C n] emission depends on three parameters (see Equations 14 
and 15): (1) the cloud column density, given by A ViC ; (2) the incident ISRF flux, xo, i n units 
of the expected flux at the Sun; and (3) the photoelectric heating efficiency, e. We have 
computed cloud models for A v c ranging from to 5 mag, and Xo ranging from 0.2 to 4.0. 
The [C n] results scale linearly with the photoelectric heating efficiency (Equation 11), so it 
was trivial to vary e. 

Figure 10 shows plots of I CII as a function of FIR for four different combinations of Xo 
and e. The curves show two types of behavior: "linear" behavior for low values of FIR; and 
"flat" behavior for high values of FIR. This is due to the difference between dust and gas 
heating mentioned above. The linear behavior occurs when A VjC is low enough that both 
dust and gas heating increase linearly with A v c . The flat behavior occurs when A v c is high 
enough for FUV heating of the gas to become negligible in cloud interiors, causing [C n] 
emission to extinguish, but not high enough for optical heating of dust grains to become 
negligible, so FIR emission is still produced. 

We have analyzed models for a range of Xo an d e values. The behavior of I CII and FIR 
in the linear regime can be approximated by: 

jlinear = F j R x (a5634 ± q 0003 ) ^ (0.267±0.001) e _ (lg) 

Errors represent formal uncertainties in fits to the grid of cloud models. The slope of the 
linear behavior (Equation 16) depends primarily on e and only weakly on xo- The slope 
should not depend at all on xo, since both I CII and the true FIR integrated intensity should 
be linearly proportional to Xo- Recall, however, that for comparison with observations we 
have approximated FIR using I 60 and I 100 (see Equation 1), i.e., we bias the FIR integral 
artificially towards short far-infrared wavelengths. As xo increases, the short wavelength 
surface brightness increases faster than linearly, since the peak of the FIR spectrum shifts to 
shorter wavelengths. Thus our FIR values should also increase with xo faster than linearly, 
while the integrated spectrum increases linearly. 
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The [C n] intensity in the flat regime can be approximated by: 

TFLAT 

— - 6 c -^ r ~ l = (149.05 ± 0.07) xo e. (17) 

10 ° erg cm z s 1 sr 1 

The flatness occurs because at high opacity all incident FUV radiation has been absorbed 
and converted to photoelectric heating (with efficiency e). The maximum possible [C n] 
intensity ought to equal the total integrated FUV intensity in the ISRF, multiplied by e: 



Xoe / J ISRF , A rfA = 1.44x 10" 4 Xoe. (18) 

</0.0912/im 



Equations 18 and 17 are indistinguishable to within 4%. We superimpose gray line segments 
on the curves drawn in Figure 10, representing the linear and flat portions of the curves, 
computed using the approximations in Equations 16 and 17. 

We define the transition value of FIR to be the value when I£™ EAR = Icn AT , i.e., 

f TP T m TRANSITION 

in ( ™ } 2 , i = (264.55 ± 0.19) x ^-^\ (19) 
10~ 6 erg cm" 2 s- 1 sr- 1 v 7 v ' 

The FIR value at which a transition between the linear and flat regimes occurs depends 
wholly on xo- We can use Equation 1 to derive a 100/zm surface brightness at the linear to 
flat transition for typical CNM gas (I 60 = 0.16I 100 ): 

TTRANSITION 

= (14.86 ± 0.01) X ( L267±0 - 001 ). (20) 

MJysr- 1 v 1 



4.4. Statistical Comparison of Data with Observations 

We can constrain e and xo by comparing our model for the [C n] and FIR emission 
of translucent clouds to the ISO observations of high-latitude clouds. First we use the 
simple equations that we derived for the linear, flat, and transition zones. Recall that the 
observations can be fit by a straight line with slope 2.5 x 10~ 2 (Equation 3). This implies 
that many of the positions observed are in the linear regime. Substitution of the slope into 
Equation 16 gives Xo °' 267 e = 0.044. The binned data, however (Figure 3), show a statistically 
significant departure from this slope at around FIR = 3 x 10~ 4 erg cm" 2 s _1 sr -1 . This 
appears to be the transition between the linear and flat regimes. Equation 19 gives Xo ~ 1-1, 
thus e ss 0.045. 

A more sophisticated way to determine e and xo f° r translucent clouds is to compare 
a grid of models to the data, and find the best fit. Recall that cloud column density, or 
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central visual extinction, is the independent variable, whereas FIR and I CII are predicted 
by the model. Given a grid of models, we can "work backwards" from the observed FIR 
surface brightness to determine A v c as a function of e and xo- Effectively FIR becomes the 
independent variable. The resulting A v c numbers yield a set of predicted [C n] integrated 
intensities, I™™ 10 ™ that we can compare with the data. We evaluate the "goodness of fit" 
as a function of e and Xo'- 

/-[PREDICTED T \ 2 

r2 / X _ ^CII L CU,i) / 9 -,n 

^ Vt, XO) — ^ < G PREDICTED p 1^.2 ' \^ L I 
^ V CII / CII 5 2 

where the sum is over all data points. This is the standard x 2 parameter — we label it L 2 here 
so as not to confuse the reader with the ISRF intensity, Xo- Here ,jPR EDICTED [ s computed the 
same way as I^f 10 ™ , using the FIR observed uncertainty range, FIRi ± a FIRti . 

Figure 11 shows contours of constant L 2 in the (e,xo) plane, representing 68.3%, 95.4%, 
and 99.7% confidence in the fit of our model to the data. These are equivalent to the 1, 2, 
and 3a confidence intervals for Gaussian statistics. The minimum L 2 point, indicated by a 
cross in Figure 11, gives the most probable values, e = 0.043 ^'[J^ an d Xo = 1-6. We derive a 
lower (— 3<T ) limit of xo > 1, but cannot derive an upper limit because the Figure 11 contours 
of constant L 2 continue indefinitely in the xo dimension. Apparently a complete transition 
to flat behavior, which determines the value of xo? does not occur in our HLC dataset. 

In terms of e and Xo, there is only marginal difference between the "warm" and "cold" 
HLC positions defined in §3. Separate L 2 analyses for the warm and cold data points yield 
e WARM = 0.053 and Xo.warm ^1-0 (3a lower limit); and e COL D = 0.045 ± 0.004 and 
Xo,cold ^ 0.8 (3a lower limit). Within the la uncertainties, the e values are identical. The 
lower limits on Xo are consistent with each other. In addition, both the warm and cold 3a 
confidence contours overlap each other nearly exactly. We conclude that translucent cloud 
positions with warm and cold IRAS 60/100 colors have nearly the same values of e and xo- 

Figure 12 displays the best fit model for I CII versus FIR, superimposed on the binned 
data from Figure 3. The observed decrease in [C n] emissivity at high values of FIR (Figure 
3) is reproduced by the model, signifying that HLCs with FIR > 3 x 10~ 4 erg cm -2 s _1 sr _1 
absorb most incident FUV photons. 

5. Discussion 

In this paper we have studied the [C n] and far-infrared (FIR) emission towards high- 
latitude translucent molecular clouds (HLCs). These observational signatures of the pro- 
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cessing of the interstellar radiation field (ISRF) by the clouds help to elucidate the role of 
dust grains in heating the gas. 

5.1. The [C n] emissivity and extinction 

The chief conclusion to be drawn from the new ISO data is that, for FIR < 3 x 
10~ 4 erg cm -2 s^ 1 sr" 1 , I CII and FIR are linearly correlated on ~ V scales (Figure 2), con- 
firming the 7° COBE FIRAS result, I CII oc FIR 095 , for the Galactic plane (Bennett et al. 
1994). The confirmation is not just qualitative: the proportionality constant is the same 
for both COBE and ISO. At low extinction the FIR integrated intensity is proportional to 
the H i column density, so the [C n] emissivity is constant for all low-extinction interstellar 
regions. 

The [C n] emissivity decreases for high extinction lines of sight due to the attenuation 
of FUV photons that cause photoelectric heating, an effect that was unnoticed in the low 
resolution COBE data. When we bin the data (Figure 3) they show a statistically significant 
decrease in slope for FIR > 3 x 10~ 4 erg cm~ 2 s -1 sr _1 . We do not have enough high 
extinction sources in our sample to be sure that I CII flattens completely, as predicted by 
our model. At the minimum, there is a link between the decreasing emissivity of [C n] and 
the transition from mostly "warm" to mostly "cold" values of I 60 /I 100 , because most warm 
positions have low FIR intensity and most cold positions have high FIR intensity. 

The COBE data do not show a decrease in [C n] emissivity for high values of FIR, 
probably because translucent HLCs are diluted inside the 7° COBE beam. The mean CO 
emission size of HLCs is ~ 1° (Magnani, Hartmann, & Speck 1996), but detectable CO 
emission from HLCs covers less than 3% of the sky (Magnani, et al. 2000). Thus on 7° 
scales, ~ 97% of [C n] emission is from purely atomic material, so the COBE survey yields 
only linear behavior for I CII versus FIR. The 70" ISO beam, on the other hand, is capable of 
resolving translucent HLCs where the column-averaged [C n] emissivity is much lower than 
in diffuse H i. 

The 60/im/lOO/im color is often used to infer either the ISRF flux (e.g., Dale et al. 
2001) or the grain size distribution (e.g., Laureijs, Clark, & Prusti 1991), but our work 
indicates that for translucent clouds variations in I 60 /I 100 can be accounted for entirely by 
variations in the extinction. Thicker clouds appear colder, even if the grain size distribution 
and the radiation field stay constant. Fitting our model to the I CII and FIR data yields the 
same values of the photoelectric heating efficiency, e, and the ISRF flux, xo, f° r sources with 
both warm and cold 60/100 colors. 
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5.2. The ISRF Flux and the Photoelectric Heating Efficiency 

The actual best fit values of the ISRF flux and the photoelectric heating efficiency are 
noteworthy. We found that the most probable bolometric ISRF flux incident on HLCs is 
Xo = 1-6, where Xo = 1 defines the ISRF spectrum of Mezger, Mathis, & Panagia (1982) 
and Mathis, Mezger, & Panagia (1983). We also derive xo ~ 1-6 independently for cloud 
MBM-12 using the observed I 60 /Iioo behavior as a function of visual extinction (Figure 5). 
The FUV portion of a xo = 1-6 field has an integrated flux of 2.9 x 1CT 2 ergcur 2 s _1 . In 
units of the Habing (1968) flux (1.6 x 10~ 2 erg cm~ 2 s _1 ), this gives G = 1.8, close to the 
actual value (Go = 1.7) determined from space-based UV measurements (Draine 1978). 
Given that our model is so successful at determining Gq, even without observing clouds in 
the flat regime where all FUV is absorbed, it is likely that the curve in Figure 12 is a realistic 
representation of the average [C n] and FIR properties of translucent high-latitude clouds. 
This means that the decrease in [C n] emissivity that we detect is indeed the transition 
between regions with less than complete and complete FUV absorption, confirming the use 
of the term "translucent" to describe the HLCs. We predict that the [C n] emissivity should 
decrease further for HLC positions with higher extinction than those in our current sample. 

The photoelectric heating efficiency that we derive, e = 4.3%, is slightly higher than the 
value 3% derived from COBE data by Bakes & Tielens (1994). Recall that e is the fraction of 
far-ultraviolet radiation absorbed by grains that is converted to gas heating. Bakes & Tielens 
estimated e for the Galaxy by dividing the observed [C n] integrated intensity by FIR. In so 
doing, they equated I CII with photoelectric heating, as we do; but they also equated FIR with 
the absorbed FUV flux, which is only partially true for translucent material. In the cold 
neutral medium half of the ISRF flux that heats dust grains and causes FIR emission is at 
optical wavelengths. The FIR integrated surface brightness can only be used to estimate the 
absorbed FUV flux if clouds have nearby O and B stars, causing FUV photons to dominate 
the radiation field. 

Another way to obtain a crude estimate of e for the cold neutral medium (CNM) is 
to consider the limiting case of complete FUV absorption in translucent gas. The [C n] 
cooling that should result is I™™ = 1.44 x 10~ 4 Xoe erg cm -2 s~ x sr _1 (Equation 18). For 
our Figure 2 HLC sample, I™* 1.1 x 10~ 5 erg cm -2 s" 1 sr -1 , so the product % e ~ 0.076. 
Taking xo — 1-6 (see above) gives e = 4.8%, implying that the CNM is extremely efficient 
at processing FUV radiation. In fact, both the 4.8% estimate and our best fit appraisal, 
e = (4.3^'2)%, are consistent with the value 4.9% adduced as the maximum efficiency for 
neutral grains (Wolfire et al. 1995; Bakes & Tielens 1994). Fewer than 50% of the grains 
responsible for heating the CNM are expected to be ionized (Li & Draine 2001; Weingartner 
& Draine 2001), so the grains should have close to maximal efficiency 
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5.3. Limitations in the Grain Model 

5.3.1. Thermal Equilibrium Assumption for Very Small Grains 

Our results depend on the assumption that all dust grains are in thermal equilibrium. 
In particular, Equations 9 and 12 require this condition. As mentioned in §4.2.1, thermal 
equilibrium breaks down for grains smaller than a ~ 25(M, which are heated by individual 
photons. These very small grains (VSGs) get much hotter than the temperatures we derived 
by equating the average power absorbed with that emitted. The effect is most obvious 
at short wavelengths (A < 60//m). Thus our predicted thermal equilibrium 60/im surface 
brightnesses are underestimates to the actual values derived when the stochastic heating of 
very small grains is considered. According to the models of Li & Draine (2001, Figure 8), 
stochastically heated VSGs are responsible for approximately half of the 60/im flux from 
diffuse high-latitude clouds (see also Verter et al. 2000). Clearly this is an overestimate 
to the stochastic VSG contribution to the emission from translucent clouds, since the power 
absorbed (and hence emitted) by VSGs is a strong function of the FUV flux, which decreases 
inside translucent clouds. If our predictions for I 60 are underestimated by a factor < 2, then 
our FIR predictions will be underestimated by < 20% (using Equation 1 and the fact that 
I 60 ~ 0.16I 100 for average H i clouds). Adding 20% to our model FIR values decreases by 20% 
the slope of the linear I CII vs. FIR relationship for a given (e,Xo) pair. Thus the derived 
product eXo°' 267 (Equation 16) decreases by 20%. This does not change significantly our 
conclusions, since e = 0.043 * 00 2 is already uncertain by 23%. We need to take into account 
the small grain physics in more detail to assess accurately the errors due to our equilibrium 
approach. This will not change an important conclusion derived from our model: variations 
in the [C n] emissivity and the 60/100 color in translucent clouds reflect variations in grain 
heating, not changes in grain optical properties. 



5.3.2. Variations in Grain Optical Properties 

It is likely that our conclusions will change for extremely opaque clouds, whose interiors 
may be sufficiently dense and free of UV photons for smaller dust particles to coagulate into 
larger grains. Two compelling pieces of evidence for variations in the optical properties of 
grains are (1) an abrupt drop in 60/100 at A v ~ 5 (e.g., Laureijs, Clark, & Prusti 1991; 
Laureijs et al. 1995); and (2) far-infrared spectra indicating an extremely cold (12 — 
13 K) population of dust towards dense cores (e.g., Bernard et al. 1999; Juvela et al. 
2002). Neither feature is a characteristic of our model, so we caution against applying our 
conclusions directly to dark clouds. Nonetheless, the effects should be negligible in the 
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model's intended domain, i.e., translucent clouds (A v = 1 — 5). 

5.4. Details of Cloud Structure 

In our analysis we made no explicit mention of the cloud density structure. We parametrized 
our model in terms of the visual extinction, which allows us to ignore density structure along 
the line of sight. If the volume density varies along the line of sight, the size scale associ- 
ated with extinction varies, but this does not affect our calculations if all else remains the 
same, i.e., the photoelectric effect is the major heating source and [C n] emission is the 
dominant cooling source. This is probably true. The brightest lines observed in HLCs 
(I CII ~ 1 x 10 -5 erg cm" 2 s" 1 sr" 1 ) account for approximately 5% of the incident FUV flux 
(when Go = 1-7; see §5.2 and Draine 1978). In other words, the FUV heating is converted to 
[C n] cooling with very high efficiency This probably requires that the [C n] excitation con- 
ditions are always thermal, i.e., either the density is above the critical density, n crit = A u i/^ ui , 
where j u i is the collisional excitation coefficient; or the temperature is above hv u ijk = 91 K. 
The average CNM temperature is ~ 80 K, but it ranges from 20-200 K (Kulkarni & Heiles 
1987), so the H i gas is probably hot enough for [C n] to be thermalized. For cooler gas 
(T ~ 20 K), the [C n] critical density is high. For excitation by H 2 , n crit m 9000cm -3 ; for 
excitation by H i, n crit ps 6700 cm" 3 . The CO-emitting molecular regions of HLCs could have 
such high densities (Ingalls et al. 2000), but it is not clear that the [C n]-emitting regions 
are so dense. It is probably the case that gas which is too rarefied to be supercritical either is 
hot enough for [C n] to be thermalized or does not contribute much to the overall extinction. 
We conclude that the detailed structure along the line of sight therefore has little effect on 
our analysis. 

Structural variations in the plane of the sky are more difficult to account for, for two rea- 
sons. First, clouds with structure in three dimensions are more permeable to the interstellar 
radiation field, especially when one takes into account scattering of FUV photons by dust 
grains (Spaans 1996). Second, the emission observed towards clouds is always the result of 
averaging the real emission across a telescope beam. Comparing two tracers observed with 
the same beamsize towards the same position minimizes the effects of structured emission. 
Unfortunately, the IRAS data we use have a resolution of 4', whereas the ISO [C n] mea- 
surements have only 71" resolution — a factor of 11 in beam area! We probably have enough 
data points that structural variations averaged out in statistics (compare Figure 2 to Figure 
3). Nevertheless, to improve the robustness of our conclusions requires higher resolution FIR 
observations, perhaps made with the Space Infrared Telescope Facility. 
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6. Conclusions 

A linear relationship is observed between the FIR and [C n] emission for high Galactic 
latitude molecular clouds observed with ISO that is indistinguishable from the prediction 
based on high-latitude COBE data, implying that the [C n] emissivity is constant for all 
low-extinction gas. At high extinction the [C n] emissivity begins to decrease due to the 
attenuation of the FUV portion of the interstellar radiation field. In contrast to other work, 
we find that differences in the 60/xm/100/mi color can be accounted for solely by extinction 
variations. Sources with both warm and cold colors seem to be exposed to the same radiation 
field, equal to the mean field near the Sun; and have the same photoelectric heating efficiency, 
close to the value for neutral grains. The transition from sources with warm to those with 
cold 60/100 colors coincides approximately with the transition from constant to decreasing 
[C n] emissivity. Such regions where the FUV spectrum is attenuated and softened, but is 
still important for heating the gas, define the translucent regime. 

We are indebted to Laurent Cambresy for providing us with an extinction map of MBM- 
12, produced using his adaptive-grid star count algorithm. This paper benefited from the 
penetrating comments and helpful suggestions of George Helou, Francois Boulanger, and 
Jean-Philipe Bernard. We thank the anonymous referee, whose comments and suggestions 
helped to improve the final manuscript. This work is based in part on the Ph.D. thesis of 
J. I. We thank Dan Clemens, James Jackson, Harlan Spence, and Antony Stark for critical 
reading of an early version of the manuscript. This research was supported in part by ISO 
under JPL contract #961513. The research described in this paper was carried out in part 
by the Jet Propulsion Laboratory, California Institute of Technology, and was sponsored by 
the National Aeronautics and Space Administration. 

A. Deriving Bulk Properties of the Grain Distribution 

We assume that grains are spherical with radius a, and come in two kinds, graphite 
(i = 1) and "astronomical silicate" (i — 2). The fundamental optical properties of interstellar 
grains — absorption efficiency, Q a bs,\( a ^)] scattering efficiency, Q sca ,x{0','i)] and mean cosine 
of scattering angle, g\(a, i) — have been computed by Draine & Lee (1984) and Laor & Draine 
(1993). The albedo can be constructed from these basic properties, oo = (Q S ca/Qabs)/[l + 

(Qsca/ Qabs)] ■ 

Following Mathis, Rumpl, & Nordsieck (1977), we assume a power law distribution of 
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gram sizes: 



drij = C;a 



-3.5 



da, a min < a < a 



(Al) 



'max t 



where drii is the number of grains per H atom of type i with radii between a and a + da, and 
the upper and lower size limits are a max ~ 0.25yum and a min ~ 0.005/iin, respectively. The 
coefficients C were determined by Draine & Lee (1984) to be C\ = 10~ 2516 cm 2,5 (H atom) -1 
for graphite and C<i = 10 -25 ' 11 cm 2 ' 5 (H atom) -1 for silicate. The values of the Cj coefficients 
were derived by fitting the theoretical extinction curve (see Equation A2 below) to the 
observed average extinction curve for the Galaxy (Savage & Mathis 1979). 

Bulk observables can be computed by integrating a particular property over the distri- 
bution of grain sizes, a and materials, i. For example, to obtain the wavelength-dependent 
extinction, A\, one integrates the extinction cross section, a extt \(a, i) = 7ia 2 Q extt \(a, i) [where 
Qext,x(a,i) = Qabs,x(a,i) + Q sca ,x(a,i)], over the grain distribution: 



The extinction curve is usually expressed as A\/A v , which requires that we multiply Equation 
A2 by the standard relationship between column density and visual extinction N H ( cm -2 ) = 
1.87 x 10 21 A v (mag) (Bohlin, Savage, & Drake 1978). 
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Fig. 1. — High Galactic latitude [C n] emission is correlated linearly with FIR emission 
on the large scale. We plot here the integrated intensity of [C n] ( 2 P 3/2 — > 2 Pi/ 2 ) emission 
from COBE FIRAS Line Emission Maps, I CII , as a function of the integrated FIR surface 
brightness, derived from COBE DIRBE ZSMA 60 and 100/xm maps resampled to the 7° 
FIRAS grid. The 5523 data points with \b\ > 5° can be fit by a straight line with slope 
(2.54 ± 0.03) x 10~ 2 . We superimpose a gray wedge on the data that represents the ±3<7 
range of expected {FIR, I 

cnj pairs. 
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Fig. 2. — The HLCs are average sources of [C n] and FIR emission. Data points show the 
ISO integrated intensity of [C n] versus FIR from IRAS observations for 101 positions in 9 
high-latitude clouds. The gray wedge labeled "COBE" represents the ±3<7 range of expected 
(FIR, I en) pairs based on fits to the all sky COBE data for \b\ > 5° (see Figure 1). The 
solid line is a weighted least squares fit to our HLC data, I CII = 2.5 x 10~ 2 FIR. 
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Fig. 3. — The [C n] emissivity of HLCs decreases as FIR increases. Dashed line seg- 
ments display weighted linear fits to the HLC data (Figure 2), separated into bins of width 
10~ 4 erg cm -2 s" 1 sr^ 1 along the FIR-axis. Circles with error bars denote the weighted mean 
values of I CII and FIR in these bins. The gray wedge represents the range of COBE data 
for average Milky Way sources, as described in the caption to Figure 2. The solid line is the 
fit to our HLC data (Figure 2). 
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Fig. 4. — Surface brightness at 100/iin versus that at 60//m for our HLC sample. A dashed 
line indicates the mean behavior derived using COBE DIRBE observations of Galactic gas 
associated with H i, I 60 = 0.16 1 100 (Dwek et al. 1997). This line is not a fit to the data. 
We identify two populations based on the 60/100 color: a "warm" population for which 
ho ^ 0.16I 1()O ; and a "cold" population for which I 60 < 0.16I 1()() . 
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Fig. 5. — The surface brightness ratio, I 60 /Iioo decreases as a function of visual extinction 
for cloud MBM-12. Image data for over 2,000 pixels are averaged in bins of width AA V = 
0.71 mag. Note that a small adjustment in 60/100 represents a large change in A v . Models 
of the thermal emission from dust grains in clouds immersed in the interstellar radiation 
field (see §4.3.3) are represented as gray curves. The contours are labeled with xo, the total 
ISRF flux measured in units of the local value. 
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MX (jam) 

Fig. 6. — The interstellar radiation field (ISRF) and the attenuated spectrum for four model 
clouds. The top curve shows the ISRF spectrum as given by Mezger, Mathis, & Panagia 
(1982) and Mathis, Mezger, & Panagia (1983), plus the 2.7 K cosmic background. This 
curve defines xo — 1- The lower curves plot the wavelength- weighted mean intensity of the 
radiation field at the centers of model clouds with central visual extinction, A v c , equal to 
0.5, 1.0, 2.0, and 5.0 magnitudes. A dashed vertical line at A = 0.2066/im, or photon energies 
of 6eV, denotes the optical/far-ultraviolet (FUV) transition. Photoelectric heating occurs 
to the right of this line. 
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Fig. 7. — The heating of an interstellar cloud with A v c = 1.0 mag. The dashed line (right 
scale) shows the equilibrium temperature of graphite grains of size a = 0.1/xm, plotted as a 
function of A v . The solid line (left scale) shows the gas heating function in the same cloud, 
assuming the photoelectric efficiency is 4.3%. We only show the profiles for half of the cloud, 
since the profiles are symmetric about A v = A v c . 
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Fig. 8. — The far-infrared (FIR) surface brightness spectra for thermal equilibrium dust emis- 
sion from model clouds. Shown here are the spectra for clouds with xo — 1 an d A v c = 0.5, 
1.0, 2.0, and 5.0 magnitudes. Vertical arrows indicate the peak of each curve, demonstrating 
that more optically thick clouds appear "colder." Thick gray curves show modified black- 
body emission for clouds with the same values of A v c , but with constant "dust temperature," 
Td =19, 18, 16.9, and 15.1 K, respectively. Dust emissivity is assumed to vary as A~ 2 . 
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Fig. 9. — The dust and gas cooling towards translucent clouds with xo — 1 and e = 4.3%. The 
two lower curves (right scale) represent emission in the 60/im and 100/im IRAS bands from 
thermal dust grains, "observed" toward a set of model clouds with different column density, 
given by A vc . The I 60 values have been multiplied by 5 for clarity. The two upper curves 
(left scale) represent the "FIR" integrated intensity, derived using I 60 and I 100 (Equation 1); 
and the [C n] integrated intensity, multiplied by 50 for clarity. 
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Fig. 10. — Observables I CII and FIR for model translucent clouds. Each curve shows the 
behavior as we vary A V)C , for four combinations of the surface radiation field, xo, and the 
photoelectric heating efficiency, e. On the 1 C11 -FIR plane, A v c increases from left to right 
and from bottom to top. The thick solid curve shows the best fit HLC values (see Figure 11), 
Xo = 1-6 and e = 4.3%. The thin solid curve shows the effect of decreasing the ISRF intensity 
to xo = 0.7, keeping e fixed at 4.3%. Dashed curves show the effect of changing the heating 
efficiency from e = 2% to 7%, holding Xo steady at 1.6. Gray line segments represent the 
linear portions of the curves and the flat portion of the Xo = 0.7 curve, calculated according 
to the approximations in the text (Equations 16 and 17). 
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Fig. 11. — Statistical comparison of observed [C n] integrated intensity and FIR surface 
brightness with our model, for a grid of e and Xo values. Contours (smallest to largest area) 
show the 1,2, and 3a confidence intervals in the fit of the model to the data shown in Figure 
2. The cross represents the best match, e = 0.043 and xo = 1-6. 
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Fig. 12. — Binned HLC measurements of I CII versus FIR (Figure 3), compared with model 
results. The gray wedge represents the range of COBE data, as described in the caption to 
Figure 2. The dashed curve represents the best fit model, e = 0.043 and Xo — 1-6. A solid 
line shows the fit to our HLC data (Figure 2). 
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Table 1. ISO LWS Observations of High-Latitude 

Clouds 



a 

bourcc Name 


(°) 


b 

(°) 


I 


FOV b 

'X') 




PA c 
\ ) 


Detected / 
Observed^ 


MBM-37 


6 


.007 


36 


.750 


21 


.0 


X 


1 


.2 


90 


7/7 


(L183) 


5 


.778 


36 


.987 


21 


.0 


X 


1 


.2 


90 


7/7 




6 


.235 


36 


512 


21 


.0 


X 


1 


.2 


90 


7/7 




7. 


172 


37 


040 


1 


.2 


X 


1 


,2 





1/1 


HD210121 


56 


.876 


-44. 


461 


1 


.2 


X 


1 


,2 





1/1 




56 


.107 


-44. 


188 


1 


.2 


X 


1 


.2 





1/1 




56 


.436 


-44 


.147 


1 


.2 


X 


1 


.2 





1/1 




56 


.511 


-44 


019 


21 


.9 


X 


1 


.2 


55 


8/9 


MBM-28 


141 


.755 


39. 


.047 


1 


.2 


X 


1 


.2 





1/1 


(Ursa Major) 


140 


.784 


38 


.501 


1 


.2 


X 


1 


,2 





0/1 




141 


.306 


39 


430 


1 


.2 


X 


1 


,2 





1/1 




142 


.002 


38. 


.700 


1 


.2 


X 


1 


,2 





1/1 




142 


.070 


38 


.199 


1 


.2 


X 


1 


,2 





1/1 


MBM-12 


159 


.251 


-34 


.483 


1 


.2 


X 


1 


.2 





1/1 


(L1457) 


159 


.151 


-33 


817 


1 


.2 


X 


1 


.2 





1/1 




159 


.301 


-33. 


501 


1 


.2 


X 


1 


.2 





1/1 




159 


.301 


-34 


471 


30 


.0 


X 


6 





135 


20/20 e 




159 


.401 


-34. 


450 


1 


.2 


X 


1 


.2 





1/1 
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.601 


-34. 


.501 


1 


.2 


X 


1 


.2 





2/2 


G225.3-66.3 


225 


.274 


-66 


.280 


1 


.2 


X 


1 


.2 





1/1 




225 


.312 


-65 


.413 


1 


.2 


X 


1 


.2 





0/1 


G295.3-36.2 


295 


.311 


-36. 


.141 


21 


.0 


X 


1 


.2 





6/7 




295 


.284 


-36 


123 


21 


.0 


X 


1 


,2 





7/7 


G300. 1-16.6 


300 


.055 


-16 


.620 


1 


.2 


X 


1 


,2 





1/1 


(Chamaeleon) 


299 


.674 


-16. 


318 


1 


.2 


X 


1 


.2 





1/1 




299 


.914 


-16. 


799 


1 


.2 


X 


1 


.2 





1/1 




300 


.250 


-16 


955 


1 


.2 


X 


1 


.2 





1/1 




300 


.263 


-16. 


767 


1 


.2 


X 


1 


.2 





1/1 


Stark 4 


310 


.142 


-44. 


475 


18 


.0 


X 


1 


.2 


90 


6/9 




309 


.578 


-45 
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1 


.2 


X 


1 


,2 
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309 


.747 
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.579 


18 


.0 


X 


1 


.2 


90 
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309 


.967 
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522 


1 


.2 


X 


1 


.2 





1/1 


MBM-33 


359 


.150 


36 
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1 


.2 


X 


1 


.2 





1/1 


(L1780) 
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.897 


36 
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1 


.2 


X 


1 
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358 


.954 


36. 


832 


1 


.2 


X 


1 


.2 





1/1 



Nomenclature for objects with |b| > 25° is that of the Magnani et al. 
(1996) catalog, with alternate names in parenthesis. Separate positions 
within the same source are listed separately. Separate observations of the 
same position, including raster maps, are listed together. Nominal source 
"center" position, if one has been designated, is listed first. 

b Ficld of view of raster map. Maps were typically made on a sparsely 
sampled rectangular grid. 

c Position angle of raster map in degrees counterclockwise from celestial 
North. 

d Number of positions detected in [C II] emission/Number of positions 
observed. 

"Previously published observations (Timmerman et al. 1998). 



